Magnetization Plateaus in a Solvable 3-Leg Spin Ladder 
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We present a solvable ladder model which displays magnetization plateaus at fractional values of 
the total magnetization. Plateau signatures are also shown to exist along special lines. The model 
has isotropic Heisenberg interactions with additional many-body terms. The phase diagram can be 
calculated exactly for all values of the rung coupling and the magnetic field. We also derive the 
anomalous behaviour of the susceptibility near the plateau boundaries. There is good agreement with 
the phase diagram obtained recently for the pure Heisenberg ladders by numerical and perturbative 
techniques. 
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One of the surprising aspects of low-dimensional quan- 
tum systems with long range finite interactions is the 
occurence of fractional magnetization plateaus. Theoret- 
ical [0-|5) considerations have revealed that such plateaus 
originate from interactions beyond those of nearest neigh- 
bours. Spin ladders provide natural examples of systems 
with non trivial magnetization plateaus since they can 
be reformulated as spin chains with longer range interac- 
tions. They have the advantage that they are accessible 
experimentally and spin ladders in a field have attracted 
quite some attention over the past few years 0. Frac- 
tional values of the total magnetization were measured 
only very recently in ID compounds and in a 2D 
system [g. 

Groundstate phase diagrams and magnetization plat- 
eaus have been calculated for a special family of ladder 
models for which partially exact results could be obtained 
p0| . We consider an integrable 3-leg spin ladder, or spin 
tube for periodic boundary conditions [ pdp2] |. Its phase 
diagram can be calculated for all values of the rung cou- 
pling and the magnetic field from the Bethe Ansatz so- 
lution. The spins along each leg and each rung have 
an isotropic Heisenberg interaction, with the introduc- 
tion of many-body terms to retain integrability. The 
model is a generalisation of Wang's 2-leg ladder model 
|p3| . The overall phase diagram of the 2- and 3-leg lad- 
ders compare well with those obtained recently for the 
pure Heisenberg ladders, by DMRG, bosonization, series 
expansions and mappings to effective Hamiltonians (see, 
e.g., Refs. Illll). 

The Hamiltonian of our model is 

H = E + E ^ runs + E ^ fiold ' C 1 ) 

i—l i—l i—l 

where, 

i=i 



H r s = E yi (*p ■ - 1) , (3) 

tff eld = -/*5>*)f- (4) 
i=l 

The operators (cr x )\ l \ and (er z )f^ act as the cor- 

responding Pauli matrices on the (i, l)th factor in the 
Hilbert space. 

As shown in |Q , the Hamiltonian ([!]) is integrable for 
h = 0. The addition of the magnetic field term how- 
ever, does not destroy integrability since [H?f, Hf cM + 
^ficidj _ q T n f n ow i n g we set J\ = Ji = J and 
J j, = J', so that we can go from the isotropic tube, 
J' = J, to the ladder, J' = 0. 

For h = it is convenient to change to the basis where 
the square and the z-component of the total spin of a 
given triangle, S — erW + <x( 2 ) + cr^ are diagonal [p"2| . 
It follows that the eight states on a given triangle fall 
into a spin- 1 quadruplet and two spin-i doublets. We 
will denote these states by |2s z ) g for the quadruplet and 
|2s 2 )di, (i = 1,2) for each of the doublets. Switching on 
the magnetic field breaks this symmetry further due to 
the Zeeman splitting. The energies of the rung and field 
Hamiltonians are given by 

H mns + ff fiold = diag{-3 J - h, -3 J + h, 

-2 J' - J-h, ~2J' -J + h, -3h, -h, h, 3h}, (5) 

on the states 

{|+k, l-k, l+k, \-)d 2 , 1 + 3)„ \+) q , \-) q , I - 3) J. 

Since S 2 and S z commute with H, the total spin and its 
z-component are good quantum numbers, as in Ref . |l0[ . 

H can be diagonalized using the Bethe Ansatz. It 
is important to note that the Hamiltonian (^) does not 
change under the change of basis given above. Further- 
more, ^ is invariant under any choice of reference state 
(or pseudo- vacuum) |f2) and any assignment of Bethe 
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Ansatz pseudo particles. For each choice however, one 
has to re-interpret this assignment. The rung and field 
Hamiltonians do alter with the choice of ordering, but 
the change is just a rearrangement of their eigenvalues 
along the diagonal. We use this property to our advan- 
tage by doing calculations with that choice of ordering 
for which the Bethe Ansatz reference state is closest to 
the true groundstate of the system. The eigenenergies of 
Ef=i H lf+i are S iven b y 



Mi 



1 



(0) 



where the number \j satisfy the well known Bethe 
Ansatz equations p7|| , 



A« + 



Mi x (l) x (l) 



AL 



A (1) 



k (2) 



and for r = 2, . . . , 7 with M s = 0, 



: ' r A (r) 



n 



x (r) 
M r _i ,(r) 

11 ,(r) 
fc=l A j 



Al r) +i 

_ \(r-l) 



M r4 



i(r+l) 



A 



(r-l) 



in 



fe=l *j 



A 



(r+1) 



(8) 



where ,7 = 1,..., M r . 

In the following we will restrict ourselves to the quad- 
rant J, h > and consider the magnetization, which is 
defined by 



M 



i=l {=1 

where n is the number of legs. 



(0 



(9) 



(M) = 1. For J — h > 2 there is another massive phase 
where the groundstate consists of singlets on each rung. 
Here evidently (M) = 0. In between these two phases lies 
a massless phase where the magnetization varies contin- 
uously. On approaching the lines h = J ± 2 from within 
the massless phase, the susceptibility shows the familiar 
square root singularity Jl9[ . 

On the line h = J the spin up triplet and the singlet 
state are degenerate. Therefore, on this line the magne- 
tization (M) = h. At the point J = h = log 2 the other 
excitations become massless and a completely massless 
phase is entered. This phase actually covers a finite re- 
gion around the origin which seems to be common to 
this type of solvable ladder model. The other point that 
can be calculated exactly marking its phase boundary 
is J = 2, h = 0. The line h = J can be regarded as 
the onset of a plateau boundary. Although there is no 
singularity in the magnetic susceptibility in the present 
case, it will appear as soon as the plateau opens. It is 
expected that the opening of this plateau is governed by 
anisotropy 




FIG. 1. The solvable 2-leg ladder phase diagram. The 
dashed line h = J divides phase space into two regions. The 
bold lines given by h = J ± 2 are phase boundaries. 



THE 2-LEG LADDER 

In this section we briefly review and expand the results 
of Wang jl3| before treating the 3-leg case. For the 2-leg 
case the rung states fall into a singlet and a triplet. The 
phase diagram in the J, h > quadrant is determined by 
the competition between the singlet state and the spin up 
state of the triplet. The difference between their respec- 
tive energies changes sign at h = J. This line therefore 
divides phase space into two regions. In each of these re- 
gions a convenient choice of Bethe Ansatz reference state 
and pseudo particles may be made. The value of the gaps 
can then easily be calculated and it follows that there is a 
massive phase for h — J > 2 where the groundstate is the 
simple product of the spin up triplet state. In this phase 



THE TUBE: J' = J 

Due to the extra symmetry in this case, the two dou- 
blets are degenerate for the rung Hamiltonian and only 
three eigenvalues of the Hamiltonian H\ nnK + H\ f eld are 
competing, —(3 J + h), —(3 J — h) and — 2th. Their dif- 
ferences change sign at the lines, 2h = 3 J and 4h = 3 J, 
dividing the phase space in three regions. In each of these 
regions we can make a convenient choice of ordering the 
states to facilitate our calculations. 

i) h > | J. We choose our states to be ordered in 
increasing energy with respect to h xung + h RcM , i.e. {| + 
3)g> |+)dn |+)da> • • •}• The energy, up to an irrelevant 
constant, is given by 
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+ 3JM 3 + (2/i - 3J)M 4 + 3 JA/ 6 + 2/iM 7 . (10) 

From this expression it follows straightforwardly that for 
2h — 3 J > 4 we have Mi = for the groundstate energy. 
In this phase, the pseudo vacuum is the true groundstate 
and the total magnetization (M) — 1 and all excitations 
are gapped. Below the line h = § J + 2, a finite den- 
sity of \+}d! and |+)d 2 states appear in the groundstate. 
This phase is massless and the mean magnetization varies 
continuously. As the magnetization reaches its maximum 
value, the susceptibility diverges as, 

X~(|J + 2-/ l )-i (11) 

Upon reaching the line h = |j, the three states have 
become degenerate in energy while all other excitations 
remain massive for large enough J. On this line the mag- 
netization is (M) = | |ll|. This line can be regarded as 
a plateau boundary in a similar way as the line h — J for 
the 2-leg ladder. The plateau will open upon introduc- 
tion of anisotropy. 

The first excitations to become massless on this line 
are \+) q , \—)di and | — )d 2 - This happens at the point 
3J = Tr/V^-log3. 

H) j J < h < I J. Here we choose our state to be 
ordered as {!+)&, \+)d 2 , \ + 3) 9 , . . .}. The energy in this 
case takes the form, 

Mi 

E = - V — m + (3 J - 2h)M 2 + {Ah - 3 J)M 3 

+ (3J -2h)M 5 + 2h(M 6 + M 7 ). (12) 

Following the same reasoning as above, we find that for 
3 J — 2h large enough M 2 — for the groundstate. In this 
region the magnetization is (M) = 4 and every magnetic 
excitation is massive. Because of the degeneracy of the 
spin up states of both doublets, there are massless ex- 
citations because the groundstate is essentially that of 
a spin-i chain. The | + 3) 9 excitation becomes mass- 
less at 3 J — 2h = 2 log 2 where the massless phase is 
entered. Here (M) may vary continuously. The suscepti- 
bility shows the square root singularity upon approaching 
the plateau from within this phase. 

Hi) h < |J. Here, the ordering for the lowest energy 
rung states is given by {]+)&, | +)<*,, \~)di, |-)d 2 > • • •}• 
The energy in this case takes the form 

Mi 

E = - V — + 2hM 2 

+ (3 J - Ah)Mi + 2h(M 5 + Al 6 + M 7 ). (13) 

Above the line h = log 2 there is the plateau phase with 
(M) = |. Below this line, the excitations \—)d t and |— )d 2 



become massless and the magnetization is allowed to vary 
continuously. On the line h = the two doublets are 
degenerate and the magnetization (M) = 0. As before, 
we may regard the line h = as a plateau boundary. 
Our results for the tube are summarised in Fig. 0. 
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FIG. 2. The solvable 3-leg tube phase diagram. The dashed 
lines 2h = 3J and 4h — 3J divide phase space into three 
regions. The bold lines h — |j + 2, h — | J — log 2 and 
h — log 2 are phase boundaries. 

THE LADDER: J' = 

For the ladder the degeneracies are lifted and the mass- 
less regions display more structure. Phase space is now 
divided into six regions due to level crossings, of which 
the boundaries are given by h = -|j for k = 1,2,3,4,6. 
The only boundaries where level crossings concerning the 
groundstate occur are those with k — 3 and k — 6, i.e. 
the same as those for the tube. A similar analysis as for 
the tube gives the following results. 

i) h > | J. As in the case for the tube, for strong mag- 
netic fields the ladder is completely magnetized. Below 
the line h = | J + 2, some of the \+)d 1 states appear 
and the magnetization may vary continuously. Again, 
the susceptibility diverges with a square root singularity 
as this plateau is approached. On the line 2h — 3 J, the 
states | + 3) q and are degenerate and the magne- 
tization thus is given by (M) = |. As before this line 
is interpreted as a plateau boundary. The line does not 
extend to the origin, but at some finite value of J = log 2 
other magnetic excitations become massless. 

ii) jJ < h < ^J. For 3 J — 2h large enough only the 
state appears in the groundstate. In this region the 
magnetization is (M) = i and every excitation is mas- 
sive. As opposed to the tube, the two doublets are not 
degenerate and the groundstate is completely polarized, 
as in Ref. [|o|. The | + 3) 9 excitation becomes massless 
at 3J — 2h = 2 where the massless phase of the previous 
paragraph is entered. Here (M) may very continuously 
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but the phase is still polarized. Also here the square root 
singularity shows up when approaching this plateau. 

On the line Ah = 3J, the (M) = § plateau extends 
into the massless phase. On this line the states | + 3) g 
and \—)di are degenerate so that they may combine to a 
net nonmagnetic excitation. It ends at some finite value 
of J which cannot be calculated analytically. 

Hi) h < | J. Here, the ordering for the lowest energy 
rung states is given by |+)d 2 , |-)d 2 , • • •}• 

Above the line h = 2 there is the plateau phase with 
(M) = i. Below this line, the excitation \—)d! becomes 
massless and the magnetization is allowed to vary contin- 
uously. For J not too small, the phase remains polarized. 
The line h = 0, where the magnetization (M) = 0, may 
again be regarded as a plateau boundary. The second 
doublet enters the groundstate at J = log 2 [[l2| and the 
polarized excitations become massless. As before, the 
susceptibility shows the square root singularity upon ap- 
proaching this line from above and we may regard the 
line h = as a plateau boundary. 

Our results for the tube are summarised in Fig. ^. 




FIG. 3. The solvable 3-leg ladder phase diagram. The 
dashed lines 2h = 3J and 4/i = 3J divide phase space into 



three regions. The bold lines h - 
boundaries. 



} J ± 2 and h — 2 are phase 



CONCLUSION 

In summary we have considered a simple solvable 
model that displays fractional magnetization plateaus. 
Such plateaus have been studied theoretically only by 
perturbative or numerical techniques. The phase dia- 
gram and in particular the location of the plateaus and 
the expected square root singularity have been calculated 
exactly. The phase diagrams in Figs. [l||| are to be com- 
pared with those of the pure Heisenberg ladders given 
in Ref. The overall agreement is seen to be excel- 
lent and our results should provide a useful benchmark 
for further studies. Consideration of the complete eigen- 



spectrum will allow the exact study of thermodynamic 
effects. More intricate phase diagrams can be obtained 
similarly for solvable n-leg ladder and tube models. 

In future work we hope to be able to include anisotropy 
in order to tune the width of the plateaus, see, e.g., the 
discussion in Ref. with respect to the XXZ chain. 
For example, in this way the lines with magnetization 
(M) = 0, (M) — | and (M) = | may disappear or may 
broaden to reveal a finite gap. The fact that the models 
in this Letter map onto an isotropic low energy effective 
Hamiltonian in a strong coupling approach E3] is consis- 
tent with this picture. 

This work has been supported by the Australian Re- 
search Council. 



[1] M. Oshikawa, M. Yamanaka, I. Affleck, Phys. Rev. Lett. 

78, 1984 (1997). 

[2] K. Totsuka, Phys. Lett. A 228, 103 (1997); Phys. Rev. 

B 57, 3454 (1998). 
[3] D. C. Cabra, A. Honecker and P. Pujol, Phys. Rev. Lett. 

79, 5126 (1997); Phys. Rev. B 58, 6241 (1998). 
[4] K. Hida, J. Phys. Soc. Jpn. 63, 2359 (1994). 

[5] K. Okamoto, Solid State Comm. 98, 245 (1996). 

[6] See, e.g., G. Chaboussant, M.-H. Julien, Y. Fagot- 

Revurat, M. E. Hanson, L. P. Levy, C. Berthier, M. Hor- 

vatic, and O. Piovesana, Eur. Phys. J. B 6, 167 (1998), 

and references therein. 
[7] Y. Narumi, M. Hagiwara, R. Sato, K. Hindo, H. Nakano 

and M. Takahashi, Physica 246B-247B. 
[8] W. Shiramura, K. Takatsu, B. Kurniawan, H. Tanaka, 

H. Uekusa, Y. Ohashi, K. Takizawa, H. Mitamura and 

T. Goto, J. Phys. Soc. Jpn. 67, 1548 (1998). 
[9] H. Kageyama, K. Yoshimura, R. Stern, N. V. Mushnikov, 

K. Onizuka, M. Kato, K. Kosuge, C. P. Slichter, T. Goto 

and Y. Ueda, Phys. Rev. Lett. 82, 3 168 (1999). 

[10] A. Honecker, F. Mila and M. Troyer, |sond-mat/9910438| , 

to appear in Eur. J. Phys. B. 
[11] M. T. Batchelor and M. Maslen, J. Phys. A 32, L377 

(1999). 

[12] J. de Gier, M. T. Batchelor and M. Maslen, Icond 



mat/9911074, to appear in Phys. Rev. B. 



[13] Y. Wang, Phys. Rev. B 60, 9236 (1999). 

[14] F. Mila, Eur. Phys. J. B. 6, 201 (1998). 

[15] K. Tandon, S. Lai, S. K. Pati, S. Ramasesha, and D. Sen, 
Phys. Rev. B. 59, 396 (1999). 

[16] R. Citro, E. Orignac, N. Andrei, C. Itoi and S. Qin, J. 
Phys.: Condens. Matter 12, 3041 (2000). 

[17] B. Sutherland, Phys. Rev. B 12, 3795 (1975). 

[18] A possible plateau at this value, i.e. (S z ) = |, was dis- 
cussed in Ref. |h|. 

[19] P. W. Kasteleyn, J. Math. Phys. 4, 287 (1963); G. I. 
Dzhaparidze and A. A. Nerseyan, JETP Lett. 27, 334 
(1978); V. L. Pokrovsky and A. L. Talapov, Phys. Rev. 
Lett. 42, 65 (1979). 



4 



